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Abstract 


The dynamic loads of a helicopter 
rotor in forward flight are influenced 
significantly by the geometric pitch 
angles between the structural axes of the 
hub and blade sections and the olane of 
rotation. 

The analytical study presented in- 
cludes elastic coupling between inplane 
and out-of -plane deflections as a function 
of geometric pitch between the plane of 
rotation and the principal axes of inertia 
of each blade. In addition to a mean col- 
lective-pitch angle, the pitch of each 
blade is increased and decreased at a one- 
per-rev frequency to evaluate the dynamic 
coupling effects of cyclic feathering mo- 
tions. The difference in pitch between 
opposed blades gives periodical coupling 
terms that vary at frequencies of one- and 
two-per-rev. Thus, an external aerody- 
namic force at n-per-rev gives forced res- 
ponses at n, n±l, and n±2 per rev. 

The numerical evaluation is based on 
a transient analysis using lumped masses 
and elastic substructure techniques. A 
comparison of cases with and without cyclic 
feathering motion shows the effect on com- 
puted dynamic rotor loads. The magnitude 
of the effect depends on the radial loca- 
tion of the pitch change bearings. 

Introduction 

For a St if f -in-plane rotor system, 
the blade chordwise stiffness may be 20 
to 50 times greater than the blade beam- 
wise stiffness. The elastic structure 
tends to bend in the direction of least 
stiffness, resulting in dynamic coupling 
between out-of -plane and inplane motions 
as a function of the geometric pitch 
angles due to collective pitch, built-in 
twist, forced cyclic feathering motions of 
the torsionally-rigid structure, and elas- 
tic deformation of the blade and control 
system in the torsional mode. 

Typical cruise conditions for a mod- 
ern helicopter require collective pitch 
angles of 14 to 16 degrees at the root, 
depending on the amount of built-in twist. 
Cyclic feathering motions of 6 to 7 de- 
grees are required to balance the one-per- 
rev aerodynamic flapping moments. In cur- 
rent design practice, elastic torsional 


deflections of the blade and control sys- 
tem of a stif f-in-plane rotor are generally 
less than one degree. The largest part of 
the angular motion in the blade-torsion de- 
gree of freedom, therefore, is the forced 
feathering motion due to cyclic pitch. 

Periodic variations of the inplane/ 
out-of-plane elastic coupling terms are 
caused when the geometric pitch angle of 
each blade is increased and decreased at a 
frequency of one cycle per rotor revolution. 
When one blade is at high pitch and the op- 
posed blade is at low pitch, an asymmetri- 
cal physical condition exists with respect 
to a reference system oriented either to 
the mast axis or to the plane of rotation. 
Oae-half revolution later, the reference 
blade is at low pitch and the opposed blade 
is at high pitch. Thus, periodic dynamic 
coupling occurs at the principal frequency 
of one-per-rev with respect to a rotating 
coordinate system. The coupling terms are 
nonlinear functions of blade pitch; hence, 
these terms also have 2-per-rev content. 

Both the steady and periodic coupling 
terms have been treated in an analytical 
study of the effects of one-per-rev cyclic 
feathering motions on dynamic rotor loads. 
Equations have been derived and programmed 
for a digital computer solution of the 
transient response of an elastic two-bladed 
rotor. 

The rotor is modeled by elastic sub- 
structure elements and lumped masses, for 
which the accelerations and velocities are 
integrated over small time increments to 
determine time histories of deflections, 
inertia loads, bending moments, etc. The 
time-variant analysis includes the capa- 
bility to calculate rotor instabilities. 

The present computer program has been 
tested for this capability, but further 
discussion of instabilities is beyond the 
intent of the paper. 

Dynamic rotor loads have been calcu- 
lated for a parametric series of rotors, 
where the coupled natural frequencies were 
tuned over the range of contemporary de- 
sign practice for teetering rotors, a de- 
scription of the analysis and a sununary of 
computed results is presented. 
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Objective 


A primary consideration in the design 
of a helicopter rotor is to minimize os- 
cillatory bending loads, or at least to 
reduce the loads to a Level that will en- 
sure satisfactory fatigue life. During 
early stages of design, the principal 
method of evaluating the dynamics of a 
proposed rotor is to calculate its coupled 
rotating natural frequencies. If required, 
design changes are made to achieve suffi- 
cient separation between the natural fre- 
quencies and harmonics of the rotor oper- 
ating speed. 

Current practice at Bell Helicopter 
Company is to require a separation of 0.3 
per rev for all flight combinations of 
rotor speed and collective pitch. One 
purpose of the present analytical develop- 
ment is to determine whether the separa- 
tion rule may be relaxed due to beneficial 
effects of cyclic feathering motions on 
rotor dynamic response. 

Collective and Cyclic Modes 

The calculation of natural frequen- 
cies for semi-rigid rotors uses a coor- 
dinate system that is based on the plane 
of rotation. The orientation of the cen- 
trifugal force field, the angular motion 
allowed by the flapping hinge(s), and the 
constraints of opposing blades lead to the 
segregation of natural frequencies into 
collective modes, cyclic modes, and (for 
four-bladed rotors) scissor or reaction- 
less modes. This procedure allows the use 
of continuous -beam theory for a single 
blade, where the centerline boundary con- 
straints are imposed from conditions of 
symmetry or asymmetry to match deflections, 
slopes, shears, and moments for the other 
blades . 

The centerline boundary conditions for 
the collective mode (Figure 1) are: 

- zero vertical (out-of -plane ) slope 
change , 

- vertical deflection constrained by 
mast tension/compression, 

- inplane slope constrained by mast 
torsion, and 

- zero inplane translation. 

The centerline boundary conditions for 
the cyclic mode (Figure 2) are: 

- vertical slope change unrestrained 
(except with flapping springs), 

- zero vertical deflection, 

- zero inplane slope change, and 


- inplane deflection constrained by 
mast shear. 



FIGURE 1. TYPICAL COLLECTIVE MODE 

FREQUENCIES AND MODE SHAPES. 



QUENCIES AND MODE SHAPES. 
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For the reactionless modes, the cen- 
terline boundary conditions are: 

- zero slope change and zero trans- 
lation in both the inplane and 

vertical directions. 

Uncoupled frequencies are determined 
by setting the geometric pitch angle of 
each elastic element to zero. The un- 
coupled frequencies are shown in Figures 
1 and 2 by the labeled curves. Note that 
the frequencies of the vertical (out-of- 
plane) modes are highly dependent on rotor 
speed, and that the frequencies of the in- 
plane modes are only slightly dependent 
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Coupled natural frequencies are shown 
as small circles in the figures. Typical 
collective modes have very small frequency 
shifts as a function of collective pitch. 
However, the cyclic modes (Figure 2) 
couple significantly with collective pitch. 
Note that the inplane frequency decreases 
and the vertical frequencies increase with 
collective pitch. The method of deter- 
mining these coupling effects is given in 
Reference 1. 

By using only one blade plus appro- 
priate boundary conditions, this method of 
calculating rotor natural frequencies is 
based on one explicit assumption, i.e., 
all other blades are at the same geometric 
pitch angle as the reference blade. If 
the blades are at different pitch angles, 
then the conditions of symmetry or asym- 
metry are not present- The inclusion of 
cyclic feathering motion, therefore, re- 
quires that the analysis treat separately 
each blade of the rotor and provide a 
means of matching the centerline slopes 
and deflections. 

Elastic Substructured Rotor Analysis 

A digital computer program has been 
developed to study the effects of cyclic 
feathering motions on dynamic rotor loads. 
The Bell Helicopter computer program is 
identified by the mnemonic ESRA for Elas- 
tic Substructured Rotor Analysis. 

The analysis is a transient solution 
of elastic rotor blade motions, where the 
coupling terms for each blade are treated 
separately and the necessary constraints 
are imposed on each blade to insure slope 
and deflection continuity at the rotor 
centerline. Each blade is considered as 
being divided into a discrete number of 
segments, with uniform weight and stiff- 
ness properties over the length of a seg- 
ment. The geometric pitch angle of each 
segment is a function of rotor azimuth 
position and input values of fore/aft and 
lateral cyclic pitch. To represent the 
hub structure that is inboard of the pitch- 


change bearings, the inboard elastic ele- 
ment may be specified as an uncoupled 
element (geometric pitch equals zero). 

All forces are applied at the ends of 
the elastic elements. Slope and deflection 
changes over the length of a segment are 
based on linear moment distributions versus 
span. For compatibility, shear over the 
segment length must be constant, which re- 
quires concentrated forces for both inter- 
nal and external forces. In its simplest 
form, the analysis follows a lumped-mass 
approach. All of the important rotor dy- 
namic characteristics may be retained with 
this method, however, by using the con- 
cepts of equivalent structural segments. 

The dynamic response equations are 
solved by a step-by-step iterative method in 
order to include transient conditions. If 
the initial deflections and velocities are 
specified (spanwise distributions for each 
blade), then the internal bending moment 
distributions arc found with respect to the 
rotating reference system. Internal shear 
distributions are obtained from the moment 
distributions, and summed with applied air- 
load forces and inertial components of the 
centrifugal force field to determine span- 
wise distributions of accelerations. 

The first estimates of deflection and 
velocity changes are calculated for con- 
stant acceleration during the integration 
time step. Then bending moments, shears, 
and accelerations are calculated for the 
end of the time step. Subsequent deflec- 
tion and velocity estimates are based on 
accelerations changing linearly with time, 
and the iterations continue until a pres- 
cribed error limit is satisfied for the 
entire set of accelerations, or until a 
limit is reached on the number of itera- 
tions . 

In recognition of the problems inher- 
ent with this type of numerical integration, 
the initial development of the ESRA com- 
puter program has been limited to a quali- 
tative study of cyclic feathering effects. 
The current program represents each blade 
with four elements, each with beamwise and 
chordwise bending elasticity. Only the 
forced rigid -body motion is allowed in the 
blade-pitch degree of freedom, i.e., elas- 
tic blade torsion is not considered. 

The current computer program is lim- 
ited to two-bladed rotors, and the torsion- 
al impedance of the drive system is assumed 
to be zero. In practice, the rotor senses 
a two-per-rev torque from the mast that is 
proportional to the drive-system impedance 
times the Hooke ’ s- joint angular oscillation, 
which is a function of rotor flapping. The 
rotor is the predominant inertia component 
of the drive system, and a good approxi- 
mation for two-bladed rotors is to assume 
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that the true axis of rotation remains per- 
pendicular to the tip-path plane even when 
the tip path plane is not perpendicular to 
the mast. Thus, Coriolis accelerations 
equal to the product of coning times flap- 
ping are not appropriate in a two-bladed 
rotor analysis. 

Bending deflections of the elastic 
elements are linearized; therefore, Cori- 
olis accelerations from radial foreshort- 
ening are excluded also. Vertical and in- 
plane translational motions of the rotor 
center are not included in the current 
version of the program. 

Referring to the description above, 
the formulation of the analysis allows the 
removal of these limiting assumptions. 

For instance, nonlinear bending deflections 
and Coriolis accelerations may be included 
by a direct addition to the inertial forces 
acting on each mass. Translation of the 
rotor centerline, additional blades, con- 
trol ^system flexibility, elastic blade 
torsion, and nonlinear hub and control 
kinematics also may be added within the 
existing computational method. 

With the limitation of four elastic 
elements for each blade, plus provisions 
for slope and deflection continuity at the 
rotor centerline, the current ESRA pro- 
gram allows 15 distinct vibration modes 
for the rotor : 

- 3 rigid -body modes (flapping, mast 
torsion, blade pitch) 

- 6 coupled elastic collective modes 
(3 vertical, 3 inplane) 

- 6 coupled elastic cyclic modes (3 
vertical , 3 inplane ) 

In attempts to predict rotor loads for 
two-bladed rotors, emphasis is placed on 
response components at least up to the 
third harmonic of rotor speed. Three-per- 
rev airloads excite the cyclic mode that 
derives from the first elastic asymmetric 
mode in the out-of-plane direction. Four 
elastic elements for each blade should 
provide a very satisfactory dynamic repre- 
sentation for this frequency range. At a 
frequency of f ive-per-rev , the second elas- 
tic mode would be excited and computed 
loads may be marginally valid. Current 
design practice is to minimize higher fre- 
quency loads by proper tuning of the rotor 
natural frequencies, as discussed earlier. 

Numerical Evaluation 

A parametric computer study was ac- 
complished to resolve a basic question: 

With respect to the natural fre- 
quency of the first coupled vertical 


elastic cyclic mode at or near 
3 per rev, how much does cyclic 
feathering motion affect 3-per- 
rev dynamic rotor loads? 

Selection of Rotor Dynamic Characteristics 

Corresponding to a Huey main rotor, 
the computer study was based on a 48-foot 
diameter 2-bladed semi-rigid rotor, oper- 
ating at 300 RPM. Two basic design ap- 
proaches were selected as end points for 
the evaluation . 

1. A constant blade weight distri- 
bution of 1.20 Ib/in. with no dynamic 
tuning weights was picked to simulate the 
early production Huey rotors. Uniform 
beamwise and chordwise stiffness values 
were determined to locate the two lowest 
coupled cyclic mode frequencies at 1.40/rev 
(inplane) and 2.60/rev (vertical) for a 
collective pitch of 14.75 degrees. The fan 
plot of cyclic mode natural frequencies for 
this rotor is shown in Figure 3. 



FIGURE 3. CYCLIC MODE, TUNED 
BELOW 3/REV. 


2. A very recent rotor development at 
Bell (the Model 645 rotor) was simulated by 
a configuration with a constant blade weight 
distribution of 1.00 Ib/in plus a dynamic 
tuning weight of 100 pounds located at the 
blade tip. Uniform beamwise and chordwise 
stiffness values were determined to locate 
the coupled cyclic mode frequencies at 
1.40/rev (inplane) and 3.40/rev (vertical), 
as shown in Figure 4, again for 14.75 de- 
grees of collective pitch. 
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FIGURE 4. CYCLIC MODE, TUNED 
ABOVE 3/REV. 

Between the two basic configurations, 
a series of intermediate rotor parameters 
was established by stepping the uniform 
blade weight from 1.20 down to I. 00 by 
increments of 0.025 Ib/in., while increas- 
ing the tip weight from 0. to 100. by in- 
crements of 12.5 pounds. Beamwise and 
chordwise stiffnesses were varied to hold 
the coupled inplane frequency at 1.40/rev 
while tuning tfxe coupled vertical fre- 
quencies from 2.60/rev to 3.40/rev in in- 
crements of 0.10/rev. Thus, to compensate 
for the program restriction of no hub mo- 
tion, proper placement of the coupled in- 
plane frequency was maintained by the 
selection of rotor stiffness. This ap- 
proach affects the spanwise distribution 
of inplane bending moments, but is entire- 
ly adequate for a qualitative evaluation. 

All of the above frequencies were 
tuned with the first segment uncoupled 
(hub structure to .25 radius), which maxi- 
mized the coupling of the vertical mode 
near 3/rev and minimized the coupling of 
the inplane mode near 1/rev. 

Additional input data was taken 
directly from the Bell Helicopter Rotor- 
craft Flight Simulation, program C81-68 
(References 2, 3), for a Model 309 King- 
Cobra flying at 150 knots. Data used in 
the present computer evaluation included 
a collective pitch setting at the root of 
14.75 degrees, a total cyclic pitch of 
6.30 degrees, and the spanwise distributed 
airloads up to and including the third 
harmonic components. 


The study results presented below are 
based, therefore, on full-scale parameters 
that are realistic with regard to current 
helicopter design practice. Although di- 
rect correlation with measured loads is not 
possible because of the simplifying assump- 
tions, it may be noted that the magnitude 
of calculated bending moments is well with- 
in the expected range. 

Computed Results 

The forced response was computed for 
the series of nine parametric rotor con- 
figurations, where the inplane coupled fre- 
quency was held at 1.40/rev and the verti- 
cal coupled frequency was varied from 2.60/ 
rev to 3.40/rev. The dynamic rotor loads 
for each configuration were calculated 
twice, once with cyclic feathering and 
once without cyclic feathering. 

Figure 5 shows the 2 /rev vertical 
bending moment at the rotor centerline as 
a function of natural frequency of the ver- 
tical elastic cyclic mode. The 1/rev vari- 
ation in structural coupling due to cyclic 
pitch, and the 3/rev applied airloads pro- 
duce a 2 /rev component of bending moment. 
This additional component peaks and changes 
sign as the vertical mode is tuned through 
3/rev. For the two-bladed rotor. I/rev and 
3/rev vertical bending moments at the cen- 
terline are negligible. 
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Inplane bending moments at the rotor 
centerline are shown in Figure 6. The 
large peak in the overall oscillatory mo- 
ment occurs as the coupled vertical mode 
is tuned through resonance at 3/rev. Note 
that the coupling associated with cyclic 
feathering increases the L/rev response by 
about 5 percent for the vertical frequency 
tuned to 2.6/rev. In other respects, the 
effect of cyclic feathering appears to be 
minimal. 


Beamwise moments and chordwise mo- 
ments at mid span are shown in Figures 7 
and 8, respectively. Two-per-rev moment 


WITHOUT CYCLIC FEATHERING 

WITH CYCLIC FEATHERING 



COUPLED VERTICAL NAT. FREQ., PER REV 
FIGURE 6. INPLANE MOMENT AT CENTERLINE 


components are not shown in the figures 
because of their small magnitudes.^ The 
significance of the cyclic feathering ef- 
fects at midspan is consistent with that 
indicated in earlier figures for the rotor 
centerline . 


WITHOUT CYCLIC FEATHERING 
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As discussed in a previous section, 
the aeroelastic effect of blade bending 
velocity was excluded from this study by 
basing the response calculations on a pre- 
scribed set of airloads. No inference is 
intended regarding the magnitude of aero- 
dynamic damping that may be associated with 
elastic bending velocities. Conversely, 
the procedure was selected so that the 
time-variant stride tural couplings could be 
studied in an analytical environment that 
does not include other sources of damping. 

The computed responses appear as un- 
damped resonances centered at 3/rev, from 
which it follows that cyclic feathering 
motions do not provide any significant 
amount of equivalent damping to suppress 
3/rev dynamic loads. Regarding the verti- 
cal cyclic mode near 3/rev, in particular, 
the effect of cyclic feathering motion 
does not provide relief for the design rule 
that requires 0.3/rev separation of coupled 
frequencies from excitation harmonics of 
rotor speed. 

The results presented above are all 
based on rotor structural simulations with 
the inboard 25 percent radius treated as 
non-feathering hub structure. This option 
of the program was selected to maximize 
the coupling (as a function of collective 
pitch) of the vertical cyclic mode near 
3 per rev. As noted, the largest change 
in rotor loads due to the inclusion of 
cyclic feathering motions was a 5 percent 
increase in inplane bending moments at the 
rotor centerline. 

The pitch-change or feathering bear- 
ings of production two-bladed main rotors 
are located typically at about 10 percent 
radius. In this respect at least, the 
above results are based on a dynamic model 
that is not representative of actual de- 
sign practice. 

To evaluate the importance of the 
radial location of the bearings, another 
set of rotor loads was computed for a 
case where the entire radius is in the 
feathering system. 

A constant blade weight distribution 
of 1.20 Ib/in with no dynamic tuning 
weights was selected, as before, to simu- 
late the early production Huey rotors. 

The structural properties of the rotor 
were modified to maintain a 1.40/rev 
natural frequency for the coupled inplane 
cyclic mode. For the modified parameters, 
the natural frequency of the coupled verti- 
cal cyclic mode is 2.87/rev. 

The computed results are shown in 
Figures 9 through 12 for the case in which 
the feathering bearings are located at 
zero percent radius- The bar graphs show 
first, second, and third harmonics plus 


overall levels of oscillatory bending mo- 
ments. The open bars are for the condition 
of no cyclic pitch, i.e., the geometric 
pitch of the elastic elements held fixed 
at the specified value of collective pitch. 
The closed bars are for the condition that 
the geometric pitch of the elastic struc- 
ture is a function of both collective pitch 
and cyclic pitch. 

Vertical and inplane oscillatory bend- 
ing moments at the rotor centerline are 
shown in Figures 9 and 10, respectively. 

The vertical moments are not changed 



1 2 3 OVERALL 

HARMONIC (PER REV) 


FIGURE 9. VERTICAL MOMENT AT CENTERLINE 



12 3 OVERALL 

HARMONIC (PER REV) 


FIGURE 10. INPLANE MOMENT AT CENTERLINE 
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significantly by the inclusion of cyclic 
pitch. However, the inplane centerline 
moments increase by 57 percent, with both 
the first and third harmonics contributing 
to the increase. 

Beamwise and chordwise oscillatory 
bending moments at 50 percent radius are 
shown in Figures 11 and 12. Most of the 
increase in the overall oscillatory mo- 
ments at mid-span is due to an increase 
in 3/rev response. 



FIGURE II. BEAMWISE MOMENT AT MID-SPAN 



FIGURE 12. CHORDWISE MOMENT AT MID-SPAN 


Due to cyclic feathering motions, a 
significant increase (57 percent) in rotor 
loads is indicated with the feathering^ 
bearings at zero percent radius; a mini- 
mal increase (5 percent) is indicated with 
the feathering bearings at 25 percent 
radius. This suggests that the radial lo- 
cation of the feathering bearings may have 
a controlling influence on the magnitude 
of the cyclic-feathering effect. Further 
study of this relationship is in progress. 

Gone lusions 

1. The cyclic feathering motions of 
a helicopter rotor cause time -dependent 
elastic coupling due to asymmetrical pitch 
on opposed blades. The effect of these 
motions on dynamic loads may be calculated 
by modeling the rotor with elastic sub- 
structure elements, by providing individual 
treatment of each blade, and by matching 
slopes and moments at the rotor centerline. 

2. Cyclic feathering 'motion of the 
elastic blade structure does not cause any 
significant damping effect on the 3-per-rev 
dynamic loads of a two-bladed semi-rigid 
rotor. The design rule requiring 0.3/ rev 
separation between coupled natural fre- 
quencies and aerodynamic excitation fre- 
quencies should not be relaxed on the 
expectation of beneficial effects from 
cyclic feathering. 

3. The inplane one-per-rev rotor loads 
of a st if f - in-plane rotor are affected sig- 
nificantly by cyclic feathering of the 
elastic structure. The magnitude of the 
effect is decreased as the feathering bear- 
ings are moved radially away from the rotor 
centerline . 
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